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Comment on "Magnetization Process of Single 
Molecule Magnets at Low Temperatures" 

In a recent letter Fernandez et al. |l| present Monte 
Carlo (MC) simulations of relaxing interacting dipolcs 
for various lattice structures and argue that the short- 
time behavior of the magnetization M(t) in an initially 
depolarized sample {M(0) = 0) is lattice structure de- 
pendent. They claim that in a simple cubic (SC) lattice 
at short times, AM(t) = \M(t)-M(0)\ follows the usual 
\ft law, but that in body centered (BCC) and face cen- 
tered (FCC) cubic lattices, AM(t) ~ t - 7 . 

We disagree with this conclusion, and we believe it 
arises from finite size effects in the MC simulations in 


The s/i short-time relaxation is expected when the 
dipolar interaction dominates over the hyperfine cou- 
pling, ie., when the spread of dipolar field bias Wd ^ £o, 
where £ is the range of energies over which the nuclear 
bias fluctuates Q (in a completely demagnetized sam- 
ple Wd is several times larger than the strength Ed of 
nearest- neighbour dipolar interactions). It should persist 
at least until t ~ (Wd/£o)t , where r Q is a microscopic 
single molecule tunneling relaxation time. In the very 
short time limit the relaxation is linear in t. Since these 
predictions all depended on the long-range part of the 
dipolar interaction, they should not depend on lattice 
structure. 

To show this, we perform MC simulations for BCC, 
FCC and triclinic lattices with and without periodic 
boundary conditions (P. B.C.). At t = all spins ori- 
entations were random and uncorrelated, with M(0) = 

and Si = ±z (\z\ = 1). The dipolar interaction 
V D {?a) = EDiSiSj-SiSinj^Sjnj)/^]/^, and dipole 
fields and spin configurations were updated at time in- 
tervals St by flipping individual spins with probability 

1 — exp{— St t^ 1 ^) exp(—£,/T)}. Here t^ 1 ^) is an in- 
coherent relaxation rate t^ 1 (^) = T ~ 1 e~\t\/to ( see |\| )• 
we put r = 1 and assume contact with a bath at tem- 
perature T. The total bias energy £ = ^ip — 9e^BSH z , 
including the internal dipolar contribution ^dip- 
In the case of BCC and FCC clusters with P.B.C., 

with dimensions 16 x 16 x 16 for Ed = 10£ o we, as in 
Q, found M{t) ~ t~°- 7 . However, we also found that 
in the FCC cluster (for example) the number of spins in 
resonance (|£| < £ Q ) is very small, much smaller than for 
the corresponding SC cluster. This means that in small 
FCC clusters, the MC procedure may not accumulate 
proper statistics. 

To solve this problem one can (i) increase the cluster 
size, or (ii) increase the fraction of spins in resonance, by 
reducing Ed/£,o- For a BCC lattice (2 molecules per unit 
cell) we applied the first scenario. For a FCC lattice (4 
molecules per unit cell) we applied both, taking Ed = 


2.5£ c . Since in a completely demagnetized FCC sample 
W D ~ 10£b, then W D /£o ~ 25, ie., still >> 1. 

The results (Fig. 1) are clear- for larger samples one 
finds \fi relaxation if Wd/£o » 1) after the usual initial 
linear behaviour. Only the duration of the initial tran- 
sient and the subsequent \Jt are sensitive to the crystal 
structure. 
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FIG. 1: M(t) vs Vt for M(0) = 0. (a) BCC lattice of size 
32 x 32 x 32, with P.B.C.; Diamonds - the MC result; solid line 
- the y/i slope of the curve; (b) FCC lattice, of size 32 x 32 x 32, 
and P.B.C. (solid line), and of size 48 x 48 x 48, no P.B.C. 
(diamonds); and (c) Triclinic lattice of 60 3 Fe — 8 molecules, 
with no P.B.C. (Diamonds - the MC result; solid line - the 
y/i slope of the curve). Each of the 8 spin 5/2 Fe +3 ions was 
oriented along the easy axis, with lattice positions obtained 
from the Cambridge Data base (easy and hard axis orienta- 
tions were taken from ;3j). All curves were obtained after 
averaging over a few hundred random initial configurations of 
spins. 
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